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Abstract 

We consider the issue of determining the D^-D^ width difference AFd ex- 
perimentahy. The current situation is reviewed and suggestions for further 
study are given. We propose a number of decay modes in addition to 
those studied in the recent E791, FOCUS and BELLE hfetime determination 
experiments. Then we address prospects for determining CF - CDS strong 
phase differences, hke 5ktt which appears in the CLEO study of K^tt^ 
transitions. We show how to extract 6k*it with CDS data and furthermore 
show when D Ki^tt data becomes available that 6kit can also be obtained. 
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I. INTRODUCTION 



In the Standard Model, effects of D^-D^ mixing are mucli smaller than those in the 
kaon, Bd and Bs systems. However, charm-related experiments of increasing sensitivity 
have been carried out, leading to ever-improving bounds on the dimensionless mixing pa- 
rameters xd = AMd/Fd and = ArD/2rD. Most recently, the E791, CLEO, FOCUS 
and BELLE collaborations have reported on attempts to detect mixing in the D-meson sys- 
tem. This has prompted discussion in the literature as to whether actual Z)-meson mixing 
(specifically a nonzero AFd) is being seen for the first time. Since a rigorous theoretical 
prediction for AFq is unlikely, experimental progress in this area is needed. In this paper, 
we discuss specific proposals for further work in lifetime difference measurements and in 
experimentally determining the strong phase 6 (which occurs between Cabibbo-favored and 
Cabibbo-doubly-suppressed decays) . 



II. MEASUREMENTS OF LIFETIME DIFFERENCES 

The E791 §, FOCUS § and BELLE § experiments study the time dependence for 
D^{t) K+K- (CP = +1 final state) and D^{t) Kn (CP-mixed final state) under the 
assumption that CP invariance is assumed. This is reasonable in view of both theoretical 
expectations based on Standard Model physics and also recent CLEO results (see Sect. III). 
If we adopt the convention that CV\D^) = +1-0") and introduce the CP eigenstates 

\D.) = ^[\D')±m] , (1) 

then \Di) is CP-even and \D2) is CP-odd. It follows from Eq. (|1|) that 

mt)) = ^[\D,{t)) + mt))] , (2) 

where |-Di,2) evolve in time with distinct masses and decay widths, 

\D,{t)) = e-'''''-"^^'^'\Dk) (A; = 1,2) . (3) 

If the K^K~ final state is overlapped with Eq. (|2|) only the \Di{t)) part contributes, leading 
to the exponential decay equation 

Tj,K{t)=Aj,Ke-^^' . (4) 

For the Ktt final state, we express the time evolution of as 

mt)) = u{t)\D') + f4t)m (5) 

where 

hit) = I E e-^^^*-i^^* f4t) = ^ E i-y^' e~^M^-h^^ (g) 
^ k=l ^ i=l 
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Then the above conditions 1, 2 imply 

TK-iT++K+n- (t) = Ak-k 



2A^^e-(^^+^^)*/^cosh 



(ri 



(7) 



The experimental conditions are such that the cosh term in Eq. ([^) is nearly unity. Thus the 
time dependence becomes exponential, allowing determination of (Fi + r2)/2. The E791, 
FOCUS and BELLE experiments measure the quantity ?/cp, 



- 1 



(8) 



and find 



Vcp 



(0.8 ±2.9 ±1.0)% 
(3.42 ± L39± 0.74)% 

(1.0 m 



(E791) 

(FOCUS) 

(BELLE) 



(9) 



Due to its superior sensitivity the FOCUS determination dominates, the net result being a 
positive value for yep of several per cent at about the two standard deviation level. 



A. Additional Decay Modes 

We urge that additional lifetime studies on CP eigenstates of the neutral D be carried 
out. It is essential to improve the statistical data base and to acquire a sensitivity beyond 
the current 2a level. Beyond that, there is still no experimental input on the pure CP = — 1 
lifetime. By using lifetimes obtained from pure CP = ±1 modes, one would be determining 
AF directly rather than comparing an average of CP = ±1 lifetimes with that of CP = +1. 

There are a number of opportunities for further study, each final state occurring in 

decay being a potential candidate. We shall discuss just a limited number of these in 
the following, citing disadvantages as well as advantages. An important subset of our list 
has modes which contain a pair of mesons, each of which is self-conjugate under the CP 
operation. If each member of the pair has spin-zero, the orbital angular momentum is S- 
wave and the CP value of the two-particle state is simply the product of the individual 
CP values. If one meson has spin-zero and the other has spin-one, then conservation of 
angular momentum requires the particles to be in a P-wave. In this case, the CP value of 
the two-particle state becomes minus the product of the individual CP values. 

Vertex identification is a key to a successful lifetime measurement. Starting from 
its (assumed known) production point, the will travel unobserved and ultimately decay 
into some final state particles. In the best case, all these are charged and the decay point 
becomes well determined. In the worst case, each primary decay product is neutral and, 
if unstable, decays itself into neutral particles. Then even with calorimetric information, 
attempting to fix a decay point is problematic. For reference, in the E791, FOCUS and 
BELLE experiments the detected modes [Ktt and KK) had just two particles (both charged) 
in each final state and the branching fractions were BK-n+ = (3.83 ± 0.09) ■ 10^^ and 
13k-k+ = (4.25 ± 0.16) ■ 10-^. § 
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Many decay modes contain neutral kaons in the final state. The neutral kaons will in 
turn decay as Kg or Ki^ mesons. For a lifetime determination measurement, a mode is 
superior to a Ki^ mode because: (i) the detection efficiency is rather larger than the J^l 
detection efficiency, so the statistics will be better for the former, (ii) the Ki^ decay occurs 
further from the decay vertex, so its background problem is more severe. Both of these 
considerations are inherent for any detector. However, since progress in dealing with Ki^ 
detection is anticipated the Ki^ modes should not be totally disregarded. To summarize. Kg 
detection is easier and can be done now whereas i^L is harder and may be done later, although 
not as well. Finally, we note that PDG listings give branching fractions for D — »• K^X {X 
denotes other final state particles) rather than for D K^ i^X. It will suffice below to use 
the approximations 

Tno^KsX — ^DO^Ki^X — -^DO^ROX ■ (10) 

These relations are not exact because decay into Ks or Ki^ is subject to interference between 
Cabibbo favored (CF) and Cabibbo doubly suppressed (CDS) modes. We discuss aspects 
of this interference in the next section. 

Now we turn to the list of additional possible modes, partitioned according to the CP of 
the final states and presented as CP = —1, CP = +1 and CP-mixed. 

Pure CP = -1 Modes 

1. Ks4>: Both the and (f) decay into charged final states, so this mode is an attractive 
one as regards particle detection. Since the (p — > K^K~ transition is a strong decay, 
it occurs right at the decay vertex. Also, the (f) has a narrow decay width. The 
branching fraction for this mode is acceptably large [B^Ks = (4-3 ± 0.5) ■ 10^^). 

2. K^uj: Although the branching fraction is respectable {Bkslu = (1.05 ±0.2) ■ 10~^), the uj 
decays predominantly via the three-body mode tt+tt^tt^ which renders it more difficult 
regarding identification of the decay vertex. 

3. Ksp^: In this case, the branching fraction is not unattractive [Bx^pO = (0.61 ± 0.09) ■ 
10~^), and the K^p^ final state would decay into all charged particles. However, the 
larger width of the (compared to the 0) makes detection relatively more difficult. 

4. KsTT^: This mode has a reasonably large branching fraction Sj^g^o = (1.06±0.11)-10~^). 
However, the presence of the 7r° hinders accurate vertex identification. 

5. K^T] and Ksi]': Both these modes are potentially interesting since the 
branching fractions are not highly suppressed {Bkst] = (3.5 ± 0.5) ■ 10"'^ and 
Brsv' — (8-5 ± 1-3) ■ lO^'^). The problem of vertex ID for a final state r] and rj' would 
resemble that for a final state u. 

Pure CP = +1 Modes 
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1. Ti^Ti~: This mode provides a clean CP = +1 signal but has the disadvantage of a small 
branching fraction (Bt^-t^+ = (1.52 ±0.09) ■ 10^'^), about three times less than K^K~. 
Also backgrounds could be a problem since since there are more tt^tt" combinations 
in a typical event (although in the rest frame, the two pions emerge back-to-back 
with larger momenta than any other final state). 

2. Ki^T] and Ki^r)': Although the branching fractions equal those for Kgr] and Kst]', de- 
tection of the presents difficulties, as discussed earher. 

3. 7r°0: Particle ID is more of an issue than for the Ks(f) mode as the neutral 7r° decays 
via the chargeless two-photon mode. Even though the branching fraction here is 
comparatively large {Bt^oks = (1-05 ± 0.11) • 10""^), it is not sufficient to compensate 
for the detection problem. Moreover, other decay modes containing tt^'s would be a 
source of background. 

4. Ksfo{9SO): This mode consists of a scalar-pseudoscalar pair in an S-wave, and has 
CP = +1 for Ks and CP = +1 for /o(980). Ahhough the Ks and /o(980) each decay 
into charged particles, the branching ratio is small (i3i^s/o(980) = (2-9 ± 0.8) ■ 10~^). 
Similar comments apply to the i^s/o(1370) final state and to the i^s/2(1270) (except 
that here the final state is D-wave) . 

5. 0p°: This mode will have positive CP provided the and p° are in an S-wave of D- 
wave state. Both decay strongly into charged particles, so the decay point will have 
four emergent tracks. The branching fraction is rather small (i3<^pO = (6 ± 3) ■ 10"'') . 

6. Ki^TT^: This final state has the same branching fraction as K^tt^, but an even greater 
detection problem due to the Ki^. In practical terms, vertex identification would be 
an insurmountable obstacle. 

Mixed CP = ±1 Modes 

For definiteness consider a Dahtz plot analysis for the neutral {Q — 0) three-body state 
(7r7r^)Q=o- There will be resonance bands corresponding to the quasi two-body modes 
{pK)q=o and (7rA'*)Q=o- Although neither the p nor the K* is a particularly narrow reso- 
nance, these decays are CKM dominant so the branching ratios are relatively large. Specific 
examples of mixed CP = ±1 modes are: 

1. 7r~^K*~: This quasi two-particle state has a large branching fraction B^+x*- = (5.0 ± 
0.4) • 10~^ and there are the two measureable K* decay modes K*~ — >■ n^K~ and 
K*~ — > TT'K^. The latter provides a rather clean three-body configuration, 7r'^{7r~Ks) 
where the parentheses stress the K*' parentage. 

2. p^K~: The largest branching fraction among all quasi two-body final states for 
decay occurs here, Bp+K- — (10.8 ± 0.9) • 10~^. The p+ decay proceeds through only 
the mode p"^ — > tt+tt*^. 

3. 7r°K*°: The branching fraction B^o^.o = (3.1 ± 0^4) • 10"^ is relatively large. The 
associated three-body configurations will be 7r°(7r~ir°) and the less useful 7r°(7r°ir~). 
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III. MEASUREMENTS OF WRONG-SIGN TRANSITIONS 

Another study which impacts on determining AFd is the CLEO experiment which 
studies the decay rate for D^{t) K^n~. This wrong-sign process can be produced both 
indirectly, from mixing followed by a CF decay, and directly, from CDS decay. The decay 
rate is given in the CP-invariant limit by 



r(t) = e 



-t 



Rd + V Rd y' ' t -\- Rm ■ t"^ 



(11) 



The i?D term arises from CDS decay, the -Rm term from mixing and the \^Rd term from 
interference between the two. We also have the definitions 

y' = y cos6 — X sm6 , x' = x cos6 + y sm6 . (12) 

The parameter 6 is the (strong-interaction) phase difference between between the CF and 
CDS amplitudes, 

<5^5g = 5-+-5+- , (13) 
where the phases and appear in the amplitudes 

MD"^K-n+ = \MDO^K-n+\ e'^-+ M D^^K+n- = \M DO^K+n-\ e'^+- . (14) 

Note that we sharpen the notation for 5 (5 — > 6^jf^) in Eq. (p!3| ) because we will encounter 
several analogous phases in our analysis. 

The CP-invariant rate formula of Eq. (pA| ) can be generalized to incorporate various 
sources of CP- violation (CPV), § 

i?D ^ (1 + ^d) (CDS) , 
y' -^y'{l + AM/2) (mixing) , (15) 
— > S^K^ + (j) (interference) , 

where Ad, Am and (f) parameterize the extent of CP violation. When the data is fit to 
include the effects of CP violation none is found. 

Am = 0.23til ± .01 , Ad = -0.01^;}? ± .01 , sin = 0.00 ± .60 ± .01 . (16) 

In the same fit one finds at 95% D.L. 



X 



(0 ± 1.5 ± 0.2)% and y' = (-2.5^^;^ ± 0.3)% (17) 



or equivalently 

\x'\ < 2.9% and - 5.8% < y' < 1.0% . (18) 

Given the present strength of the CLEO signals for x' and y', it is prudent to cite the 
results as bounds as in Eq. (pISl). One expects future experiments to reduce the statistical 
and systematic uncertainties. Even so, ignorance of the phase S^^^ will hamper efforts to 
compare the FOCUS/E791/BELLE results with those from CLEO. 
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A. On the Determination of y' 



Can theory alone provide the value of 5/^^''? Symmetry considerations are of only limited 
use. It is known that b^^^ vanishes in the SU(3) invariant world P,|TO|, and this result has 
been recognized |]1T| in discussing aspects of the wrong-sign transitions. Thus, calculating 
the value of b^^} necessarily involves the physics of SU(3) breaking. Unfortunately, our 
limited understanding of physics in the charm region (especially the complicating effects of 
QCD) makes it difficult to perform reliable calculations. [0 It is, perhaps, not too surprising 
to find rather different statements in the literature about depending on the underlying 



approach. In one analysis |]13[, the findings of Ref. and Ref. [|T5| are shown to imply 



rather small values for b^-^ ■, less than 15°. However, the resonance model of Ref. [|T6l has 
considerably greater SU(3) breaking and obtains values as large as b'^^ ~ 30°. The largest 
value cited for b^^} appears in Ref. which shows that accepting the central values of the 
FOCUS and CLEO experiments leads to b^^^ in the second quadrant. However, it has been 
argued [jl^ that within a reasonable range of SU(3)-breaking parameters it is not possible 
to arrive at very large values of b'^^} (45° or larger) of the type considered in Ref. 

In view of this state of affairs, it makes sense to explore what experiment can teach us. 



B. Doing without Data 

Recalling our comments in the previous section on the relative measureability of -ft's.L 
modes, we begin by assuming that only a data set not containing final state -R^l's is available. 
The inclusion of data is covered later. 

Our first conclusion concerns the CDS D —>■ Ktt decays: K^tt^, —>■ K^tc^, 

K^TT^ and —>■ K^tc'^. At present, only the first of these has been observed 
{Bk-^+ = (1.46 ±0.30) ■ 10"^). If only Ks data is used, then neither D° K^tt° nor D+ 
K^TT^ modes can be determined experimentally. This can be understood by considering CF 
and CDS transitions having a neutral kaon in the final state: 



D K^ir CF Decays D K^n CDS Decays 



Transition 


Final State 


Transition 


Final State 



















For both CDS transitions —>■ K^tt^ and D'^ K^tt~^, the will decay via the same 
Ks TT+TT^ mode as the CF transitions —>■ K^tt^ and K'^tt'^. Since the CF 

decays dominate, extracting information about CDS final states containing a from just 
Ks detection will be impossible. This negates performing a direct experimental measurement 
oibS- 

What is the situation for other possible final states like Kp or K*tt1 Clearly, the same 
no-go result will hold for the D Kp decays. This leaves only the case oi D ^ K*7r. 
Since the K* decays strongly into two different charge combinations of Kir, each D K*ii 
transition will have two final configurations. Continuing to assume that only the Ks mode 
in K^ and K^ is observed, we obtain the following list: 
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D ^ K*n CF Decays 



Transition 


FSl 


FS2 


DO ^ ^*o^o 

D+ K*^7T+ 




(i^svr-)7r+ 
(irsvr°)7rO 
(i^s7r°)7r+ 



D K*TT CDS Decays 



Transition 


FSl 


FS2 


DO K*+7T- 


(ir+7r0) 


TT 


(i^svr+)7r- 


DO ^ i^*o^o 


{K+n- 


)7r0 


(irsvr°)7rO 


D+ K*+nO 








D+ K*^Tr+ 


{K+n-] 


7r+ 





Each fCvr arising from K* decay is enclosed in parentheses and FSl, FS2 are the two 
three-body final states per D decay. Each CDS transition with a K*^ in the final state 
has a configuration (FS2) identical to that of a CF transition with a K*^ in the final state. 
However, the other configurations (FSl) each contain a charged kaon and thus distinguish 
between CF and CDS decays. 

Thus, all four D K*'k CDS decays can be utilized. In those final states containing a 
K*^ , both configurations FSl and FS2 will have a unique signature (it is, however, necessary 
to employ a Dalitz plot analysis to properly identify which 'ii'vr' composite is a product of 
K* decay). For final states with a K*^, there will be a reduction factor of 2/3 in the number 
of events since only the configuration FSl can be used, i.e. 



and 
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(19) 



Thus we are led to analyze the phenomenology of D 
and CDS sectors. 



K*iT transitions in both the CF 



1. Cabibbo Favored (K*Tr) Decays 

There are three Cabibbo favored (CF) D —>■ K*tt decays, 

D° ^ K*-n+ , D° K*\° , D+ K*%+ . (20) 

These proceed through the QCD-corrected AS = AC = ±1 weak hamiltonian, which takes 
the form p[ 

= c_7^L^^) + cX''"^ • (21) 

TY^^'' and T^^^"* transform alike under isospin, as the J3 = +1 member of an isotriplet. 
Under SU(3), however, belongs to 6©6* and Tif^^ to 15©15*. [|l|,|Ol The coefficients 

c± encode the short distance, perturbative part of the QCD corrections. At energy scale 
Mw, c± have essentially equal magnitudes. As the energy scale is lowered, the coefficient c_ 
is enhanced whereas c+ is suppressed. 

Using just the isospin property of H^^\ we express the above decay amplitudes asQ 



Equivalent formulae can be written for Kn and Kp. 
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MK*On+ 



(22) 



where the subscripts represent twice the isospin of the final state K*tt composites. Observe 
that these amphtudes obey the sextet-dominance constraints A^x*o^+ = and J^k*-tt+ = 
— v^A^^.o^o. Upon either expanding the relation = {K*^'k^\[I^,H^^'^]\D^) or utilizing 



the amplitude relations of Eq. (^2]), one arrives at the isospin sum rule pT|^4 
Of interest to us here are the phase difference and amplitude ratio, 



(23) 



cos {61 - 63) 

As 
Ai 



er 



V2_ 

4 [r/f*o^+ {3TK*-n+ + 3r^*o^o — r 



K*Ot, 



11/2 



2T 



The most recent data compilation | 



1/2 



gives for the three CF modes D K*n, Kn, Kp, 



(24) 



Mode 




A3/A, 


K*7r 




U-^J-O.Ol 


Ktt 


90.2 _g 20 


0.37 ±0.03 


Kp 


0.0° ± 44.9'' 


0.39 ±0.10 



The preceding equations can of course be used to obtain phase relations in addition to 
those in the above table, e.g. for the Kn system. 



5, = 79.5° 
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5, = 110.3° 



(25) 



and so on. Staying temporarily with the D — > Kit decays, let us compare physics of the real 
world with that of a world which is SU(3) symmetric and in which c+ = 0: 









Hypothetical world 


0.5 


0.0 


Real world 


0.551 ±0.006 


0.296 ±0.028 



For these rates, at least, the agreement between the real world and the hypothetical world is 
not unreasonable. That the hypothetical world is, in some sense, nearby the real world will 
be useful later as a guiding principle in our study of the CDS amplitudes. This comparison 
between the real world and the hypothetical SU(3) world having c+ = explains the small 
observed values of A^j Ax. In an SU(3) symmetric world, the limit c+ = would correspond 
to As/Ai = for the D Kn amphtudes. Although the precise values of c± depend on the 
renormalization scheme (involving both the choice of operator basis and of renormalization 
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scale /i), a typical numerical value is c^/c^ ^ 0.5 for the range 2.0 > /i(GeV) > 1.5. 
The short distance effects embodied in c± account for much of the suppression for A^/Ai 
observed in the above table, the rest arising from the operator matrix elements. Operator 
matrix elements play a much larger role in the kaon system (A/ =1/2 rule) where QCD 
effects are more powerful. 

2. Cabibbo Doubly Suppressed (K*Tr) Decays 

Corresponding to the three D K*tt CF decays of Eq. (|20|) are the following four 
D K*n decays in the Cabibbo doubly suppressed (CDS) sector, 

D° ^ K*+n- , D+ /r°7r+ , K*\° , D+ K*+n° . (26) 

The CDS weak hamiltonian has AS = — AC = ±1 and is written analogous to Eq. (|2T|), 

H^^) = c_7^L^°^) + cX"""'^ , (27) 

but now tyI^^^-* and 7i^^^^ behave differently under isospin, transforming respectively as 
an isosinglet and as the I3 = member of an isotriplet. Under SU(3), t-^^^^^^^ (like 7i^^^) 
transforms as a member of 6©6* and ?i+ (like ) transforms as a member of 15©15*. 
Performing isospin decompositions of the CDS decay amplitudes yields 

Mk*o^o = -Aae'~^' - 2A3e'~^' , MK*+nO = -Abe'~^' + 2^36^^^ ^ (28) 

where the CDS isospin moduli and phases are labelled with super-bars. Corresponding to 
the above four decay CDS decay amplitudes are the four physical observables Aa, Ab, A3 and 
61 — 63. Two distinct J = 1/2 moduli {Aa and Ab) occur because there are two independent 
sources of the J = 1/2 final state, the isoscalar 7-^^°^^ and the isovector 

Aa = VSA^^ + aS+^ and Ab = VSA^^ - A^^^ . (29) 
We will need to determine the phase A 

A = 61-63 (30) 

and the moduli ratios r, R, 

r = ^ and R=^ ■ (31) 

Aa Aa 

Since Aa, Ab and ^3 are moduli, we have r > and R > 0. In addition, we note that 
A3/Ab = r/R. 

Taking the absolute square of each relation in Eq. (^8]) and forming ratios gives 
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7^l 
7^3 



2 + 2r2 - 4r cos A 
1 + 4r2 + 4r cos A ' 
2 + - 4r cos A 



7^2 



2 + 2r^ - 4r cos A 
2i?2 + 2r2 + 4ri?cosA 



i?2^4^2_4^^cosA ' 
where the {TZk} are the ratios of CDS decay rates, 

^ DO^K*+tt- ^ _ ^ DO^K*+it- 



7^l 



7^2 



r 



r 



(32) 



(33) 



By ehminating r and cos A from the relations in Eq. (|52D, one obtains a cubic equation 
in the variable R. However, there is an unphysical root R = —1, leaving the solution as a 
root of the quadratic equation 

7^27^3(27^l - i)r^ + (27^27^3 + 7^l7^3 - 2nin2 - 7^l7^27^3) R + 7^l7^2 - 27^l7^3 = o . 

(34) 

It turns out that to obtain the physical solution it is necessary to choose the square root of 
the discriminant as positive, 



R 



-b + \/\? - 4ac 
2a 



(35) 



where 



h = 27^27^3 + 7^l7^3 - 27^l7^2 - 7^l7^27^3 
c = 7^l7^2 - 27^l7^3 . 

To see why, let us consider a hypothetical world with 0+ = 0. Then since Ti. 
isoscalar operator, it follows that 



(CDS) 



^3 = =^ r = and 
As a consequence, we have 

ni = 2, 7^2 = 



1 



Ah = A 



7^3 



R= 1 



(36) 



IS an 



(37) 



(38) 



from which the physical solution of Eq. (p^ ) is identified. Returning to the real world, from 
R one obtains r 



(2 - 7^l)(l + i?7^2) - 2(1 - 7^2i^^)(l + 7^l: 
2(1 - 7^2)(l + 7^l) - 2(1 - 27^l)(l + i^7^2) 



1/2 



(39) 



and lastly cos A, 



cos A 



2 - 7^l + 2r2(l - 27^l; 

MiT^O 



(40) 
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3. Determining the Phase 5 



(ch) 
K'tt 



In the relation 6^^^^ = 6 



{K^ tt) tt)^ ^y^^ qp p]-^g^gg ^Cf^ pg^j^j^Q^ determined from the 
relations in Eq. (p2D because one cannot know the individual values of 6i and 63. However, 
the first relation in Eq. (^) does allow one to solve for — ^3, 



tan 



sin {61 - 63 



cos (61-63)+ A3 /2Ai_ 
Analogously, from the CDS relation in Eq. (^) we can solve for 6^_ '^^ 

sin (61 - 63 ^ 



tan 



cos 5i -^3 -^3 Ma 



where 6^_ is the phase of the CDS K*^tx amplitude. 

Combining these two relations we find for 



r(ch) 



63 — 63 + tarT^ 



sin {61 - 63 



cos {61 - 63) + A3/2A1 



— tan 



sin (61 - 53) 



cos 5i -^3 -A3/Aa 



(41) 



(42) 



(43) 



Under the assumptions that only Ks data is used and that isospin is a valid symmtery, we 
conclude that Eq. ( ^31) will be the best one can do in a purely experimental determination 
of 6^^%. An expression for 5^^*^ itself is obtained only via dropping the contribution ^3 — ^3. 
One might argue that these phases occur in exotic channels and should be individually small. 
Unlike the case of the D — > Kit amplitudes, there is no SU(3) prediction that = 0. 



C. Including Kj^ Data 

In the previous subsection, the avoidance of Ki^ data forced us to work with D K*7r 
decays. The inclusion of i^L data allows us to return to the D — > Ktt and D Kit decays 
in the following. Each D -R's.ltt mode will receive contributions from both CP and CDS 
sectors. Writing the transition amplitudes for — > i^Ts lvt'' and Kg l^t^ in a generic 

notation, we have 



Md^k, 



V2 



McfW'''^ - |-McDs|e 



Li" 



V2 



lA^crle^'^^ + |-McDs|e 



The corresponding decay widths will each contain three terms, 

r_D^Xs7r = l^CF — \/rcF\/rcDS COS {6c¥ — ^CDs) + 2-'^CDS 



Td-^Xltt = 7^C¥ + \/rcF\/rcDS COS (5cF " ^CDs) + 7^ 



CDS 



(44) 



(45) 
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To our knowledge it has been standard in the PDG data compilation for D K^n to 
ignore all but the CF contribution by using the 'factor of two rule' in Eq. (10) to infer 



the CF D — > K^X branching fraction from that oi D ^ K^X. However, some account 
of sub-dominant terms is made in Ref. by attributing to their neglect a source of 10% 



systematic error. 

In terms of Cabibbo counting, the contributions on the right hand side of Eq. (^) go as 
1 : 9l : 9^ oi roughly 1 : 0.05 : 0.002. Taking the sum of decay rates gives 

'^D^KsTT + '^D^K'L-K = TcF + TcDS • (46) 

Since no existing facility can deliver 0.2% sensitivity, the Fcds contribution to this equation 
is negligible and one arrives at the kind of relation given earlier in Eq. (^0|) . 

There is, however, the possibility of observing the 0{9l) interference term via the asym- 
metry measurement 



^ = ^— p -2—=^ COS ((5cF - ^cDs) , (47) 

J- D^KsTT + I- D^A'ltt V J- CF 

or more specifically 



Aoo = -2W cos {6ko^o - ^Ko-^o) , A+ = -'2^\\^r-^^^ cos ((5^0^+ - (5ico^+) . 

V ^ DO^RO-kO V J- D+^K07r+ 

(48) 

These asymmetries are 0{9'^), so signals will occur at about the 5% level. Such measurements 
are difficult for existing B-factories but hopefully can be performed. 

The detection of these asymmetries is clearly intriguing because they refer directly to 
6cF — ^CDS- Although the phase differences Sj^o^o — Sko^o and ^7^0^+ — 5ico^+ in Eq. ( ^8]) 
are for neutral modes (and not the charged case 5x-7r+ — ^K+n-) it is nonetheless valuable 
information. At a rigorous level, it follows from the positivity of decay widths that a negative 
(positive) asymmetry would correspond a phase difference in the first (second) quadrant. 
Beyond that one is forced into modelling Fcds- Since this contributes as a square root, 
the effect of model dependence is somewhat softened but still may be large. We note that 
^Ro-kO ~ ^ro-kO = in the SU(3) hmit. 



We conclude this section by considering how to implement a complete data set for the 
D —>■ Ktt decays. It is understood from the preceding discussion that we organize the 

and Ki^ final states into sums and differences. There will be a total of seven D — > Ktt 
decays, of which three provide information on CF physics and four on CDS physics. Defining 
Ffc = pAl/{87im'jj) {k = 1, 3), we have for CF-related decays 

^R-n+ = Fi + y^FiFg cos {61 - 63) + ^Fg 

Tk^^o + Tk^^o = ^Fi - y'FiFg cos {61 - 63) + ^Fs + F^ + A^/fJ^cos [61 - 63^ + 4f 3 

~ ^Fi - v^cos (5i - 53) + ^r3 

rR,^+ + Tks.+ = 1^3 + 2Th + 4v/rScos (5i - 63) + 2T3 ~ ^F3 . (49) 



12 



In the latter two relations, we have made the approximation of discarding 0{6f) contribu- 
tions (in accordance with the discussion around Eq. (0))- The approximate relations are 
seen to reproduce the content of Eq. (p^. For the CDS-related decays we have 

Tk+.- = - 4^/TJ'scos {Si - 63) + 2f 3 
Tk+^o = f - 4y^fbf3 cos (El - 4) + 413 

TKiTro - '^kstto = -V2\/rifaCos (^i - 5i) - 2v^y^r\f^cos (^1 - 53) 

+y2^r^cos (^3 - ^1) + 2\/2^r^cos (^3 - 4) 

rif^,r+ - rj^g^+ = 3v^\/r^cos (^3 - + 3V2\/r^cos (^3 - 4) , (so) 

where = pA|/(87rm|,) {k = a,b,3). Each term in the first two relations is 0(^c) while 
each term in the latter two are 0{9l). As an aid to analyzing these equations, we propose a 
simplified scenario with ^3 = ^3 = and = 7r/2. The approximate equations which result 
are 

Tk+^- ~ 2Ta - 4v/fJ^cos5i + 2f 3 

~ f, - A^fJlcosh + 4f 3 

r^L^o - r^^s-" - -v^\/r\flsin5i + ^^/r^cos^i + 2V2^Jv^3 

Tk^^+ - r^s-+ - Sv^^r^cos^i + 372^1^ . (51) 

In general, one must solve numerically for the unknowns Fq, Ffe, F3 and sin(5i. Let us 
point out, however, the quahtative difference between the limiting cases 5i ~ 7r/2 or 5i ~ 
0. To study this difference, it is not enough to measure just the K^tt~ and K^tt^ final 
states; the -R'l.sti" modes are required as well. It suffices to note here for 5i ~ 7r/2 that 
{T Kj^^o -T K^^o)"^ /T K+^~ ->_Fi and {T Kj^^+ -T K^ir+f k+itO F3r3/rf„ whereas for 5x^0 
both ratios become T^. As 5i proceeds from 7r/2 to 0, the first ratio decreases but the second 
increases by almost an order of magnitude. 



IV. CONCLUSIONS 

The recent FOCUS experiment on AFd has yielded a signal at the several per cent level. 
By comparison, this experimental result is over an order-of-magnitude larger than the value 
?/cp — 0.8 ■ 10~^ obtained in a theoretical analysis |]15[ based on a sum over many decay 



modes. In this paper, we have avoided the temptation to provide a theoretical prediction 
of our own for AFd- As stated earlier, we are not aware of any analytic approach in the 
charm region for which theoretical errors/uncertainties can be controlled. We therefore feel 
that whether or not the FOCUS result holds up over time is for future experimental work 
to decide. 

At the very least, however, the E791, FOCUS, BELLE and CLEO studies serve to 
stimulate fresh thinking on a subject {D^ mixing) that has long resisted progress. Our work 
in this paper has been to suggest further experimental work which would be of value: 
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1. We have described in Sect. II both positive and negative aspects of various decays 
beyond those used in the E791, FOCUS and BELLE experiments. In particular, we 
recommend that the Ks4>, K^uj and Kq,p^ modes be given serious attention. Each of 
these hes within the CP = — 1 sector, which heretofore has only been probed indirectly 
via the mixed-CP case of the D —>■ {K^tt^ + K^n^) transition. 

2. In Sect. Ill we divided our discussion of the strong phase 6 = 6^^^ into two parts: 

Supposing that accurate data on Ki^ final states is not forthcoming, we concluded 
that it will not be possible to probe the phase Sxn experimentally, but that the Sk*^ 
decays would be accessible. Thus, we propose that branching fractions for the four 
CDS decays D° ^ K*^^- , K*^tt'^ and ^ K*+7r^ , K'^7r*+ be studied. At present, 
there is data only for the K*^7i^ transition, with a stated uncertainty of about 

44%. Although any CDS branching fraction will be very small, the availability of 
copious charm production at B-factories and hadron colliders allows for the study of 
this hidden corner of charm physics.0 

We explored the eventuality that accurate data on Ki^tt final states will also be 
gathered. In principle, the asymmetries of Eq. p8|) would provide direct examples of 
CP - CDS phase differences, but are hindered by the dependence on CDS branching 
fractions. A more ambitious program would be to collect the complete set of CDS 
Kir data displayed in Eq. (0). In principle, this would allow for a determination of 
S^+J"^ like that given in Eq. (^) for K*7r. Pinally, we note that although our approach 
in this paper has been limited to what can be learned from just decay rates, the 
study of Dalitz distributions in multibody final states offers a separate opportunity for 
attacking the '5^^^^ problem'. 
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^After this paper was completed, an announcement appeared of a new CLEO measurement, 
ycp = —(1-1 i 2.5 ± 1.4)%. This is consistent with previous results (c/Eq. @). They also 
report a first measurement of the CDS mode K'^tt~t:^ which is a start of the exploration of 

the K*TT CDS modes. 
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